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1. Introduction 



The Higgs boson [1] is a fundamental ingredient of the Standard Model (SM), but it 
has not yet been observed. 

Direct searches at LEP imply a lower limit of M H > 112.3 GeV (at 95% CL) 
[2] on the mass M# of the SM Higgs boson. Global SM fits to electroweak precision 
measurements favour a light Higgs (M# ^200 GeV) [3]. The combination of the 
preliminary Higgs boson search results of the four LEP experiments [4, 5] shows an 
excess of candidates, which may indicate the production of a SM Higgs boson with 
a mass near 115 GeV. The final analysis of the LEP data is expected soon, but it is 
unlikely that it can substantially change these results. 

After the end of the LEP physics programme, the search for the Higgs boson 
will be carried out at hadron colliders. Depending on the luminosity delivered to the 
CDF and DO detectors during the forthcoming Run II, the Tevatron experiments can 
yield evidence for a Higgs boson with Mh < 180 GeV and may be able to discover 
(at the 5(7 level) a Higgs boson with M H £ 130 GeV [6]. At the LHC, the SM Higgs 
boson can be discovered over the full mass range up to Mh ~ 1 TeV after a few 
years of running [7]. 

The dominant mechanism for SM Higgs boson production at hadron colliders is 
gluon-gluon fusion through a heavy-quark (top-quark) loop [8] . At the Tevatron, this 
production mechanism leads to about 65% of the total cross section for producing a 
Higgs boson in the mass range M H = 100-200 GeV [6]. At the LHC [9], gg fusion 
exceeds all the other production channels by a factor decreasing from 8 to 5 when 
Mh increases from 100 to 200 GeV. When Mh approaches 1 TeV, gg fusion still 
provides about 50% of the total production cross section. 

QCD radiative corrections at next-to-leading order (NLO) to gg-fusion were 
computed and found to be large [10, 11, 12]. Since approximate evaluations [13] of 
higher-order terms suggest that their effect can still be sizeable, the evaluation of the 
next-to-next-to-leading order (NNLO) corrections is highly desirable. 

In this paper, we perform a first step towards the complete NNLO calculation. 
We use the recently evaluated [14] two-loop amplitude for the process gg — > H and the 
soft-gluon factorization formulae [15, 16, 17, 18] for the bremsstrahlung subprocesses 
gg — > Hg and gg — > Hgg, Hqq, and we compute the soft and virtual contributions to 
the NNLO partonic cross section. We also discuss all-order resummation of soft-gluon 
contributions to next-to-next-to-leading logarithmic (NNLL) accuracy. 

We use the approximation M t ^> Mh, where M t is the mass of the top quark. 
The results of the NLO calculation in Ref. [12] show that this is a good numerical 
approximation [13] of the full NLO correction, provided the exact dependence on 
M H /M t is included in the leading-order (LO) term. We can thus assume that the 
limit M t 3> M H continues to be a good numerical approximation at NNLO. 

The hadronic cross section for Higgs boson production is obtained by convoluting 
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the perturbative partonic cross sections with the parton distributions of the colliding 
hadrons. Besides the partonic cross sections, the other key ingredients of the NNLO 
calculation are the NNLO parton distributions. Even though their NNLO evolution 
kernels are not fully available, some of their Mellin moments have been computed 
[19] and, from these, approximated kernels have been constructed [20]. Recently, the 
new MRST [21] sets of distributions became available' 1 ', including the (approximated) 
NNLO densities, which allows an evaluation of the hadronic cross section to (almost 
full) NNLO accuracy. 

We use our NNLO result for the partonic cross sections and the MRST parton 
distributions at NNLO to compute the Higgs boson production cross section at the 
LHC. In this paper, we do not present numerical results for Run II at the Tevatron. 
Inclusive production of Higgs boson through gluon-gluon fusion is phenomenologi- 
cally less relevant at the Tevatron: it is not regarded as a main discovery channel, 
because of the large QCD background [6]. 

The paper is organized as follows. In Sect. 2 we define the soft- virtual approx- 
imation for the cross section and present our result for the corresponding NNLO 
coefficient. In Sect. 3 we discuss soft-gluon resummation for Higgs production at 
NNLL accuracy, and we also consider the dominant contributions of collinear origin. 
In Sect. 4 we present the quantitative effect of the computed NNLO corrections for 
SM Higgs boson production at the LHC. Finally, in Sect. 5 we present our conclusions 
and we comment on Higgs boson production beyond the SM. 



2. QCD cross section at NNLO 

We consider the collision of two hadrons hi and h 2 with centre-of-mass energy y/s. 
The inclusive cross section for the production of the SM Higgs boson can be written 

as 

a(s,Mfj)=y2 dxi dx ± f a/hl (xi, fi 2 F ) f b / ha (x 2 , fi 2 F ) dz5\z — 

a,b J ° J ° V X 1 X 2 

■ a z G ab (z; a s (/4), M 2 H /fi 2 R ; M 2 H /fx 2 F ) , (2.1) 

where r = Mfj/s, and fip and hr are factorization and renormalization scales, respec- 
tively. The parton densities of the colliding hadrons are denoted by f a /h{x,H 2 F) an d 
the subscript a labels the type of massless partons (a = g, qf, q~f, with Nf different 
flavours of light quarks). We use parton densities as defined in the MS factorization 
scheme. 

1"We thank J. Stirling for providing us with the new set of distributions. 
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From Eq. (2.1) the cross section a ab for the partonic subprocess ab — > H + X at 
the centre-of-mass energy s = X1X2S = Mjj/z is 

a ab (s, M 2 H ) = J a Q M 2 H G ab {z) = a z G ab {z) , (2.2) 

where the term 1/s corresponds to the flux factor and leads to an overall z factor. 
The Born-level cross section cq and the hard coefficient function G ab arise from the 
phase-space integral of the matrix elements squared. 

The incoming partons a, b couple to the Higgs boson through heavy-quark loops 
and, therefore, o"o and G ab also depend on the masses Mq of the heavy quarks. The 
Born-level contribution o"o is [8] 



G F ^ /4M 



ct = 



where Gf = 1.16639 x 10 5 GeV 2 is the Fermi constant, and the amplitude Aq is 
given by 




Aq(x) = -x l + (l-x)f(x 



arcsin* * 



(2.4) 



. 1 + y/l-X . 

In . — iTT 

1 - y/T^X 



, X < 1 



In the following we limit ourselves to considering the case of a single heavy quark, 
the top quark, and Nf = 5 light-quark flavours. We always use M t (M t = 176 GeV) 
to denote the on-shell pole mass of the top quark. 

The coefficient function G ab in Eq. (2.1) is computable in QCD perturbation 
theory according to the expansion 

G ab (z; a s (/4), M]j/ fx%) = a 2 ^ 2 R ) £ <?£>(*; 

r,—(\ \ ' 



n=0 

ry 2 (n 2 )G {0) (z) + MD (_, . M H . M 2 H 



+ — "2— G ab Z \ -T- -2" + °i a S) , ( 2 - 5 ) 

where the (scale-independent) LO contribution is 

g2 ) («) = ^^*(1-«). (2.6) 

The NLO coefficients G^ are known. Their calculation with the exact depen- 
dence on M t was performed in Ref. [12]. In the large- M t limit (i.e. neglecting 
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corrections that vanish when M H /M t — > 0) the result is [10, 11] 

<?<!> (*; = 5(1 ~z)(j + 6C(2) + In ^ + 12 2* 

In* 11 (l-^) 3 



/4 



P s 7(z) In 



'1-^) 2 M : 



2 



-6- 



1-z 



(2.7) 



MlK; M|/^) = \ Pgq (z) In (1 j,^ + \* ~ 

£ XIX -p O z 



(2.8) 



Gg^Ml/^Ml/^^gS— G«(,;M^;M^)=0, (2.9) 

where £(n) is the Riemann zeta-function (C(2) = tt 2 /6 = 1.645 . . . , C(3) = 1.202 . . . ), 
and we have defined 



VAz) 



W(l - z) 



(2.10) 



The kernels P a b( z ) are the LO Altarelli-Parisi splitting functions for real emission, 
'1 



2 + z(l-z) 



P (2) = il±(i^>! 



(2.11) 



and, more precisely, P gg s (z) is the regular part (i.e. after subtracting the 1/(1 — z) 
soft singularity) of P gg (z). 

In Eqs. (2.7)-(2.9) we can identify three kinds of contributions: 

• Soft and virtual corrections, which involve only the gg channel and give rise 
to the 5(1 — z) and T>i terms in Eq. (2.7). These are the most singular terms 
when z — > 1. 

• Purely-collinear logarithmic contributions, which are controlled by the regular 
part of the Altarelli-Parisi splitting kernels (see Eqs. (2.7), (2.8)). The argu- 
ment of the collinear logarithm corresponds to the maximum value (g^ ma![ ~ 
(1 — z) 2 Mjj/z) of the transverse momentum of the Higgs boson. These 
contributions give the next-to-dominant singular terms when z — > 1. 



• Hard contributions, which are present in all partonic channels and lead to finite 
corrections in the limit z — > 1 . 



4 



The terms proportional to the distributions T) i (z) and 5(1 — z) can be used to 
define what we call the soft-virtual (SV) approximation. In this approximation only 
the gg channel contributes and we have 



ab 



(1)SV 




+ 



33-2iV 



6 




(2.12) 



The SV terms are certainly the dominant contributions to the cross section in the 
kinematic region near threshold (r = M\j s ~ 1). At fixed s, this means that the SV 
terms certainly dominate in the case of heavy Higgs bosons. However, these terms 
can give the dominant effect even long before the threshold region in the hadronic 
cross section is actually approached. This is a consequence of the fact that the 
partonic cross section a(s, Mjj) has to be convoluted with the parton densities, and 
the QCD evolution of the latter sizeably reduces the energy that is available in the 
partonic hard-scattering subprocess. Thus, the partonic cross section a(s, Mjj) (or 
the coefficient function G(z)) in the factorization formula (2.1) is typically evaluated 
much closer to threshold than the hadronic cross section. In other words, the parton 
densities are large at small x and are strongly suppressed at large x (typically, when 
x — > 1, f(x, jj?) ~ (1 — x) v with i] ^ 3 and r] ^ 6 for valence quarks and sea-quarks or 
gluons, respectively); after integration over them, the dominant value of the square 
of the partonic centre-of-mass energy (s) = (xiX2)s is therefore substantially smaller 
than the corresponding hadronic value s. Note, also, that this effect is enhanced, in 
gluon-dominated processes, by the stronger suppression of the gluon density at large 
x. In the case of Higgs boson production at the LHC, these features were emphasized 
in Ref. [13], where the authors pointed out that the SV approximation gives a good 
numerical approximation (see also Sect. 4) of the complete NLO corrections down to 
low values (M H ~ 100 GeV) of the Higgs boson mass. 

The NNLO coefficients G ab are not yet known. Their computation, including 
their exact dependence on M t , is certainly very difficult, since it requires the evalu- 
ation of three-loop Feynman diagrams. 

The computation is certainly more feasible in the large- M t limit, where one can 
exploit the effective-lagrangian approach introduced in Ref. [22] and developed up 
to 0(a\) in Refs. [23, 13]. Using this approach, the contribution of the heavy-quark 
loop is embodied by an effective vertex, thus reducing by one the number of loop 
integrals to be explicitly carried out. 

Within the effective-lagrangian formalism, an important step has recently been 
performed by Harlander [14], who has evaluated the two-loop amplitude for the pro- 
cess gg — > H by using dimensional regularization in d = 4 — 2e space-time dimensions. 
The two-loop amplitude has poles of the type l/e n with n — 4, 3, 2, 1. The coefficients 
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of the poles of order n = 4, 3, 2 had been predicted in Ref. [24]. The agreement 
with this prediction is a non-trivial check of Harlander's result. 

To compute the NNLO cross section, the two-loop amplitude for the process 
gg — > H has to be combined with the phase-space integrals of the squares of the one- 
loop matrix element for the process gg — > Hg and of the tree-level matrix elements for 
the processes gg — > H gg and gg — > Hqq. We have computed these matrix elements 
in the limit where the final-state partons are soft, by using the one-loop and tree- level 
factorization formulae derived in Refs. [15, 16] and Refs. [17, 18], respectively. Then, 
we have carried out the phase-space integrals by using the technique of Ref. [26] . The 
result contains e-poles and finite terms. The e-poles (including the single pole 1/e) 
exactly cancel those in the two-loop amplitude [14] , thus providing a non-trivial cross- 
check of our and Harlander's results. The remaining finite terms give the complete 
soft and virtual contributions to the NNLO cross section. 

Details of our calculation will be presented elsewhere [27]. In this paper we limit 
ourselves to presenting the final result. We obtain the following soft and virtual 



contributions to the NNLO coefficient function G 



(2). 

gg 



G gg ^ SV (z] M%/ ff R , M 2 H j f/ F ) = 5(1 — z) 

19 



11399 133^ , 9 w x2 165 w x 
1u + - C < 2 > " 20 C(2) " - C(3) 



+ £>o 



+ V 1 



+ [ Y + 3 N ln M? 

(33-2N f ) 2 .pt F 



H + N 



f 



+ ~ ttt^- 1^ ^ - 18 C(2) In 2 



1189 5 A/rtN 5 A/oX 
W-3 C(2) + 6 C(3) 
Ml 



+ 



+ 



48 

169 171 
—r- + -TT- 



fi 2 R /sp 

1 9 \ M 2 

4 2 aS)--N f + (SS-2N f) a2))ln-f 

^ + f^^(33-2, /)C (2))lnf 



H + 33C(2) + |ic(3) + iV / (^-2C(2)) + 



4 2 f J & 



3 , , , M 2 M 2 

" In " 



- - (33 - 2Nf) In , 
z jj, F 



133 _ 90 C (2) - —Nt + 36 In 2 + (33 - 2N f ) ( 2 In ^ - 3 In **f) 



+ ^2 

+ 72D 3 . 



M 2 1 

-33 + 2iV f + 108 In —f 
/4 



(2.13) 



Note that our result in Eq. (2.13) gives the complete soft contributions (all the 
terms proportional to the distributions T>i (z)) to the NNLO coefficient functions 
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G ab (z). It also gives the complete virtual contribution (the term proportional to 
8(1 — z)) to the gg channel. The expression in Eq. (2.13) is an approximation of the 
exact NNLO calculation in the sense that it differs from G^ b \z) by terms that are 
less singular when z — > 1. More precisely, in the large- ,2 limit we have (see Sect. (3)) 

G $ (*) " G ?9 SV ( Z ) = OQ^O- " *)) > ( 2 - 14 ) 
G$(z) = 0(\n*(l-z)) , (2.15) 

G$(z) oc 5(1 -z) + 0((1 - z) ln 2 (l - z)) , G$(z) = 0((1 - z) ln 2 (l - z)) . 

(2.16) 

Note also that, unlike the NLO term G^ b \z), the NNLO coefficient function 
G^ b (z) is not independent of M t in the large- M t limit. The virtual contribution 
in Eq. (2.13) contains a term, proportional to In M^/M 2 , that derives from the 
integration of the heavy-quark degrees of freedom in the effective lagrangian [23, 13]. 

Our result in Eq. (2.13) can be useful as a non-trivial check of a future complete 
calculation at NNLO. It can also be used to extend the accuracy of the soft-gluon 
resummation formalism to NNLL order (see Sect. 3). 

As previously discussed, the SV approximation turns out to be a good numer- 
ical approximation of the full NLO correction for Higgs boson production at the 
LHC. Thus, the NNLO-SV result in Eq. (2.13) can also be exploited to obtain an 
approximate numerical estimate of the complete NNLO correction (see Sect. 4). 



3. Soft-gluon resummation at NNLL accuracy 

The soft (and virtual) contributions ala^V m (z) (with m < 2n— 1) to the coefficient 
function G gg (z) can be summed to all orders in QCD perturbation theory. Using 
the soft-gluon resummation formulae that are known at present, we can check the 
coefficients of some of the soft contributions presented in Eq. (2.13). The remaining 
coefficients can then be used to extend the accuracy of the resummation formulae to 
NNLL order. Both points are discussed in this section. 

The formalism to systematically perform soft-gluon resummation for processes 
initiated by qq annihilation and gg fusion was set up in Refs. [28, 29, 30, 31]. Soft- 
gluon resummation has to be carried out in the Mellin (or iV-moment) space. The 
iV-moments Gn of the coefficient function G(z) are defined by 

G N = I dzz N ~ 1 G(z) . (3.1) 
Jo 

In iV-moment space the soft (or threshold) region z — > 1 corresponds to the limit 
iV — > oo, and the distributions T> m (z) lead to logarithmic contributions, V m (z) — > 
ln m+1 N. The singular contributions in the large- iV limit can be organized in the 



7 



following all-order resummation formula: 

G m ,n = ^«(«s(^,M|/^;J^ +0(1/N) . 

(3.2) 

The radiative factor embodies all the large contributions IniV due to soft radi- 
ation. The function C gg (as) contains all the terms that are constant in the large- N 
limit and has a perturbative expansion analogous to Eq. (2.5): 



C gg (a s (fi 2 R ),M 2 H/^R,M 2 H /fi 2 F ) 

2 \ \ n 

2 ' ' • ■ " ■ - • 



1-oo / / 2 \ \ n 

i+E (^) ^; } «/^^/^) 

n=l ^ ' 



(3.3) 



These constant terms are due to virtual contributions, and the perturbative coeffi- 

— (n) 

cients C gg are thus directly related to the coefficients of the contribution propor- 
tional to 5(1 — z) in G g n g (z). The term 0(1/ N) on the right-hand side of Eq. (3.2) 
denotes all the contributions that are suppressed by some power of 1/N (modulo 
IniV enhancement) when N — > oo. 

The radiative factor A^ for Higgs boson production has the following general 
expression [28, 29, 32]: 

A> s (/4),^lr//4;^//4) = [& N (as(»l),M 2 H /faM 2 H /v 2 F )}* 

•A$*> H (a s (^),MfA£) • (3.4) 

Each term A^ embodies the effect of soft-gluon radiation emitted collinearly to the 
initial-state partons and depends on both the factorization scheme and the factor- 
ization scale In the MS factorization scheme we have the exponentiatedresu.lt 



UT z N ~ l - 1 f K ' H da 1 "i 
dz — / -\ A a (a s (q 2 )) \ , 

o i — z q ) 

(3.5) 

where A a (as) is a perturbative function 

^s) = ^4 1} +(^) 2 4 2) + p) 3 4 3) + ^s) • (3-6) 

The factor A^ clt ' ) is independent of the factorization scale and scheme and contains 
the contribution of soft-gluon emission at large angles with respect to the direction 
of the colliding gluons. It can also be written in exponentiated form as 

dz — — D H (a s ((l - z) 2 M 2 H ))} , (3.7) 
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where the function D H (a>s) for Higgs production has the following perturbative ex- 
pansion: 

D B ( aa ) = {^) 2 DP+0(<4). (3.8) 

The coefficients and fully control soft-gluon resummation up to next- 
to- leading logarithmic (NLL) accuracy [28, 29, 32]. In the case of a generic incoming 
parton a, they are given by 

A^=C ai A^ = ±C a K, (3.9) 

where C a = Cp = 4/3 if a = q, q and C a = Ca — 3 if a — g, while the coefficient K 
is the same both for quarks and for gluons [33, 30, 34] and it is given by 

k Mts-t)->>- <»»> 

Expanding the resummation formula (3.2) up to C?(a|) and transforming the result 
back to ^-space, it is straightforward to check that we correctly obtain the soft 
contributions, V (z) and V x (z), to G ( ^ SV (z) in Eq. (2.12). By comparison with the 
virtual term in Eq. (2.12), we can also extract the coefficient C^il in Eq. (3.3): 



99 

^ ) (^/^;^/^) = y + 6C(2) + ^^ln^ . (3.11) 



Then, we can expand the resummation formula (3.2) up to 0(a^), and we can com- 
pare the result with our NNLO soft-virtual calculation in Eq. (2.13). It is straight- 
forward to check that the knowledge of A^g\A^ and C gg predicts the coefficients 
of T> 3 (z),T> 2 (z) and T>i(z) in G g 2 <J SV (z), and that the prediction fully agrees with our 
result in Eq. (2.13). 

The comparison* at O(a^) and our calculation of the V -term in Eq. (2.13) also 
allows us to extract the (so far unknown) coefficient D$ that controls soft-gluon 
resummation at NNLL order. We obtain 

dP = c\ (-§ + H c(2) + 7 C(3) ) + CaNi g _ ? c(2) ) . (3 . 12) 

Note that the corresponding NNLL coefficient for the Drell-Yan process [35] differs 
from by the simple replacement of colour factors C F — > Ca- This could have 
straightforwardly been predicted from the general structure of the soft-factorization 
formulae at 0(a|) (see Sect. 5 of Ref. [16] and the Appendix of Ref. [18]). The 
exact expression of the remaining NNLL coefficient A g 3 ^ is still unknown, but an 
approximate numerical estimate can be found in Ref. [35]. 

*We can also extract the virtual coefficient C in Eq. (3.3). 
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The integrals over z and q 2 in Eqs. (3.5) and (3.7) can be carried out to any 
required logarithmic accuracy (see Refs. [32, 35]) and used for phenomeno logical 
analyses. Quantitative studies of soft-gluon resummation effects for Higgs boson 
production are left to future investigations. 

3.1 Collinear-improved resummation 

In Ref. [13] Kramer, Laenen and Spira (KLS) exploited the resummation formalism to 
obtain approximate expressions for the NNLO corrections to Higgs boson production. 
Their resummation formula is a simplified version of Eq. (3.2) that includes only the 
first-order coefficients (the coefficients A^\C^ and the first-order coefficient (3q in 
the expression of the running coupling a§(q 2 )). Therefore, the NNLO expressions 
obtained in Ref. [13] correctly predict only the coefficients of the contributions V 3 
and T> 2 to the soft and virtual coefficient function G g 2 g SY in Eq. (2.13). 

KLS also pointed out [13] that the resummation formalism can be extended to 
include subdominant contributions in the large- z limit. These contributions are the 
terms proportional to powers of ln(l — z) that appear in G gg (z) (see, e.g., Eq. (2.7)). 
In iV-moment space, they lead to contributions of the type ln fc N, which are usually 
(and consistently) neglected within the soft and virtual approximation (i.e. in the 
limit N — > oo). 

We agree with KLS that the highest power§ of ln(l — z) at the n-th perturba- 
tive order, namely, ln 2n_1 (l — z) in G g n g \z) (or, equivalently, the term ln 2n_1 iV 
in G g n g ^ N ), can correctly and consistently be implemented in the all-order resumma- 
tion formula (3.2). The key observation [13] is that these terms have a collinear 
origin. They arise from the transverse-momentum evolution of initial-state collinear 
radiation up to the maximum value of qx permitted by kinematics. In the large- 
z limit, the maximum value is qT m & x ~ (1 — z) 2 Mfj, which is very different from 
the typical hard scale Mfj of the process. The large transverse-momentum region 
(1 — z) 2 Mfj < q\ < Mfj is thus responsible for the leading ln(l — z)-enhancement. 
The resummation formalism correctly embodies the transverse-momentum evolution 
of soft radiation up to the kinematical limit (1 — z) 2 M^ (see Eq. (3.5)). Therefore, 
the leading collinear enhancement can be taken into account by supplementing the 
integrand in Eq. (3.5) with the regular (i.e. non-soft) part of the Altarelli-Parisi 
splitting function (see Eq. (2.11)). Both for the qq annihilation (Drell-Yan process) 
and gg fusion (Higgs production) channels, we can simply perform the following 
replacement on the right-hand side of Eq. (3.5): 



z 



JV-l 



- 1 



z 



N-1 



- 1 




1-Z 



A? 




A? + z 



- 1 



z n-i A m + (i/N 2 ) 



(3.13) 



l-z 



§ As for lower powers, KLS acknowledge [13] that their result is not complete. 
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Having performed the replacement of Eq. (3.13) in A^, we can insert its ensuing 
collinear-improved expression in Eq. (3.2). The resummed expression for the N 
moments of the coefficient function G g9t N can then be expanded in powers of as in 
the large-iV limit by consistently computing and keeping all the terms of the type 
a l a s~k In 2 ™ -1 N. Transforming the result back to z-space, this procedure gives the 
soft and virtual contributions to Ggg(z) plus its leading subdominant correction (the 
contribution proportional to ln 2ra_1 (l — z)) when z — > 1. 

We name soft-virtual- collinear (SVC) approximation this improved version of the 
SV expressions in Eqs. (2.12) and (2.13). We find 

Gj^^; Af^/^;^) = Gg^; JI4K;« - 12 Ml - , (3-14) 

G^{z;]tf H /fo]tf H /M . (3.15) 

The coefficient of ln(l — z) in Eq. (3.14) correctly reproduces that obtained by the 
exact NLO expression in Eq. (2.7). The coefficient of ln 3 (l - z) in Eq. (3.15) agrees 
with that computed in Ref. [13]. 

The numerical study of Ref. [13] shows that the effect of the contribution ln(l — z) 
at NLO is not small (see also Sect. 4), in particular at low values of the Higgs boson 
mass. Therefore, in our estimate (Sect. 4) of the NNLO corrections to Higgs boson 
production at the LHC, we consider both the SV approximation in Eq. (2.13) and 
the SVC approximation in Eq. (3.15). In the gg partonic channel we thus neglect 
NNLO of the type 

G$(z)-GW™(z)=0(\n*(l-z)) , (3.16) 

Note, however, that the coefficient function of the gq channel still contains con- 
tributions proportional to ln(l — z) at NLO (see Eq. (2.8)) and to ln 3 (l — z) at NNLO 
(see Eq. (2.15)). We do not consider the latter. At low values of the Higgs boson 
mass their effect is small, because the parton density luminosity of the gq channel 
is smaller than that of the gg channel. The effect increases by increasing the Higgs 
boson mass. 

4. Numerical results at the LHC 

In this section we study the phenomenological impact of the higher-order QCD cor- 
rections on the production of the SM Higgs boson at the LHC, i.e. proton-proton 
collisions at y/s = 14 TeV. We recall that we include the exact dependence on M t in 
the Born-level cross section a (see Eq. (2.3)), while the coefficient function G ab (z) 
is evaluated in the large-M t approximation. At NLO [12, 13] this is a very good 
numerical approximation when Mh < 2M t , and it is still accurate to better than 
10% when M H & 1 TeV. 
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Unless otherwise stated, cross sections are computed using the new MRST2000 
[21] sets of parton distributions, with densities and coupling constant evaluated at 
each corresponding order, i.e. using LO distributions and 1-loop as for the LO 
cross section, and so forth. The corresponding values of ^q CD (as(M z )) are 0.132 
(0.1253), 0.22 (0.1175) and 0.187 GeV (0.1161), at 1-loop, 2-loop and 3-loop order, 
respectively. In the NNLO case we use the 'central' set of MRST2000, obtained 
from a global fit of data (deep inelastic scattering, Drell-Yan production and jet 
Et distribution) by using the approximate NNLO evolution kernels presented in 
Ref. [20]. The result we refer to as NNLO-SV (SVC) corresponds to the sum of 
the LO and exact NLO (including the qg and qq channels) contributions plus the 
SV (SVC) corrections at NNLO, given in Eq. (2.13) (Eq. (3.15)). The LO and 
NLO results obtained by using the CTEQ5 distributions [36] are very similar to 
the ones computed with the MRST2000 sets (the differences are smaller than the 
uncertainties arising, for instance, from scale dependence). Therefore, we will not 
show those results ^ . 

The comparison between different sets of parton distributions, however, cannot 
be regarded as a way to quantitatively estimate the uncertainty on the parton distri- 
butions. The theoretical and experimental errors that affect present determinations 
of the parton distributions are typically larger [38] than the differences between the 
parton distribution sets provided by different groups [21, 36, 37]. In the case of Higgs 
boson production at the LHC, the study of the CTEQ Collaboration [39] recommends 
an uncertainty of about ±10% on the corresponding gluon-gluon and quark-gluon 
parton luminosities. 

We begin the presentation of our results by showing in Fig. 1 the scale dependence 
of the cross section for the production of a Higgs boson with Mh = 115 GeV. The 
scale dependence is analysed by varying the factorization and renormalization scales 
by a factor of 4 up and down from the default value Mh- The plot on the left 
corresponds to the simultaneous variation of both scales, /ip = /ir — X^h, whereas 
the plots in the centre and on the right correspond, respectively, to the results of 
the independent variation of the factorization or renormalization scale, keeping the 
other scale fixed at the default value. 

As expected from the QCD running of a s , the cross sections typically decrease 
when fip> increases around the characteristic hard scale Mh- In the case of variations 
of hf, we observe the opposite behaviour. In fact, the cross sections are mainly 
sensitive to partons with momentum fraction x ~ 10~ 2 , and in this x-range scaling 
violation of the parton densities is (moderately) positive. As a result, the scale 
dependence is mostly driven by the renormalization scale, because the lowest-order 
contribution to the process is proportional to «§, a (relatively) high power of as- 

^Larger deviations (for instance, the NLO cross section increases by ~ 10% for Mh = 100— 200 
GeV) appear when comparing to the GRV98 distributions [37], where both the gluon distribution 
and the value of a s (M z ) are different from those of MRST2000 and CTEQ5. 
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Figure 1: Scale dependence of the Higgs production cross section for Mh = 115 GeV at 
LO, NLO, NNLO-SV and N NLO- SVC. 

Figure 1 shows that the scale dependence is reduced when higher-order correc- 
tions are included and, in the case of the factorization-scale dependence, a maximum 
appears at NNLO-SV and NNLO-SVC, showing the improved stability of the result. 
Also note that there is an increase in the scale dependence when going from NNLO- 
SV to NNLO-SVC. This is due to the fact that the dominant collinear terms included 
in the SVC approximation give a sizeable contribution and are scale-independent (see 
Eqs. (3.14) and (3.15)), so their effect cannot be compensated by scale variations. 
Similar results are obtained for higher masses, with a reduction in the scale depen- 
dence when approaching high mass values. 

The impact of higher-order corrections is usually studied by computing K- 
factors, defined as the ratio of the cross section evaluated at each corresponding 
order over the LO result. The i^-factors are shown in Fig. 2 where the bands ac- 
count for the 'theoretical uncertainty' due to the scale dependence, quantified by 
using the minimum and maximum values of the cross sections when the scales fiR 
and fi F are varied (simultaneously and independently, as in Fig. 1) in the range 
0-5 < Xi XRi Xf < 2. The LO result that normalizes the i^-factors is computed at 
the default scale Mh in all cases. 

The plot on the left-hand side of Fig. 2 shows the uncertainty at LO and compares 
the exact NLO result with the NLO-SV and NLO-SVC approximations. In the case 
of light Higgs production, the NLO-SV approximation tends to underestimate the 
exact result by about 15 to 20%, whereas the NLO-SVC approximation only slightly 
overestimates it, showing the numerical importance of the term ln(l — z) added in 
the SVC approximation. Nevertheless, all the results agree within the theoretical 
bands: this confirms the validity of the large- z approximation to estimate higher- 
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Figure 2: K-factors for Higgs production for the full NLO result and the NLO-SV, NLO- 
SVC, NNLO-SV and NNLO-SVC approximations. 

order corrections, and, in particular, allows us to assume that a similar situation 
occurs at NNLO. As expected, the agreement between the three results improves for 
larger masses. 

The right-hand side of Fig. 2 shows the SV and SVC results at NNLO. Again, 
the SVC band sits higher than the SV one, the ratio of the corresponding cross 
sections being almost the same as the one at NLO, as shown in the inset plot. The 
contribution from non- leading terms ln fe (l — z), with k < 3 (which are not under 
control within the SVC approximation), is not included, but it is expected to be 
numerically less important". 

As is well known, the customary procedure (that we also are using) of varying the 
scales to estimate the theoretical uncertainty can only give a lower limit on the 'true' 
uncertainty. This is well demonstrated by Fig. 2, which shows no overlap between 
the LO and NLO bands. However, the NLO and NNLO bands do overlap, thus 
suggesting that the perturbative expansion begins to converge from NNLO. Note 
also that the size of the NNLO bands is smaller than that of the NLO bands: the 
scale dependence at NNLO is smaller than at NLO. 

Considering the results obtained at NLO, it is reasonable to expect the full NNLO 
i^-factor to lie inside the SV and SVC bands, and most probably, closer to the SVC 
one. In particular, for a light Higgs boson (M# ^200 GeV), this expectation would 
correspond to an increase of 15 to 25% with respect to the full NLO result, i.e. a 

II We have tried to add a term ln 2 (l — z) with a coefficient as large as that of the term ln 3 (l — z), 
finding only a small (about 5%) modification. 
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factor of about 2.2 to 2.4 with respect to the LO result. Taking into account that 
the NLO result increases the LO cross section by about 90% our result anticipates a 
good convergence of the perturbative series. 
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Figure 3: Cross section for Higgs boson production at the LHC in the NNLO-SV and 
NNLO-SVC approximations. 

In Fig. 3 we present the NNLO-SV and SVC cross sections as a function of 
the Higgs mass and including the corresponding uncertainty bands computed as 
defined above. To facilitate the comparison with other calculations and more refined 
predictions, we report the values of the cross sections for the production of a Higgs 
boson with M H = 115 GeV. The NNLO-SVC band corresponds to <r = 43.51-58.56 pb 
(50.13 pb at the default scales), the NNLO-SV to a = 37.73-45.69 pb (41.66 pb at 
the default scales), whereas for the full NLO it is a = 34.14-48.48 pb (40.37 pb at 
the default scales). 

Finally we want to quantify the effect of the (approximated) NNLO parton dis- 
tributions in the gluonic channel. In Fig. 4 we study this effect for the NNLO-SV 
result at Mh = 115 GeV, by plotting the cross section as a function of the scale. 
We use different combinations of NNLO and NLO parton distributions and coupling 
constant expressions. The inset plot shows the ratio R of the different results with 
respect to the one obtained by using NNLO distributions and 3-loop as- The use 
of NNLO distributions and 3-loop as reduces the NNLO cross section by 10% with 
respect to the result that would be obtained if NLO distributions and 2 loop as were 
used. Since the values of as(Mz) from MRST2000 are very similar at 2 and 3 loops 
and the typical scale of the process is not far from Mz, the effect of going from 2 to 
3 loops as amounts to only 1/3 of the 10% change. The biggest effect comes from 
the difference in the distributions, mostly due to the decrease of the NNLO gluon 
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Figure 4: Cross section for Higgs production with Mh = 115 GeV computed using different 
NLO and NNLO parton distributions and coupling constant. 

density at small x [21]. Similar results are obtained in the SVC approximation and 
for different masses. 

5. Conclusions 

In this paper we have studied the QCD corrections to Higgs boson production through 
gluon-gluon fusion in hadronic collisions, within the framework of the large- M t ap- 
proximation. Using a recent result for the two-loop correction to the gg — > H 
amplitude [14] and the soft-factorization formulae for soft-gluon emission at 
[15, 16, 17, 18] we have evaluated the soft and virtual QCD correction to this process 
at NNLO (SV approximation). We have also considered [13] the leading ln 3 (l — z) 
contribution from the collinear region (SVC approximation). Our result for the co- 
efficient G gg in Eq. (2.13) is consistent with the present knowledge of soft-gluon 

(2) 

resummation at NLL accuracy; it also allows us to fix the NNLL coefficient D H in 
Eq. (3.8). 

We have then studied the phenomenological impact of our results at the LHC by 
using the (approximate) NNLO set MRST2000 of parton distributions [21]. We have 
shown that the exact NLO result lies in between the NLO-SV approximation and 
the NLO-SVC approximation, the latter being a better numerical approximation in 
the case of low values of the Higgs boson mass. Comparing the results in the SV and 
SVC approximations at NNLO for a light Higgs (M H ^200 GeV), we estimate that 
the NNLO correction will increase the NLO result between 15 and 25%. 
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The results presented here are a first consistent (though approximate) estimate 
of QCD corrections to Higgs boson production through gg fusion at NNLO and will 
eventually be a stringent check of a future full NNLO calculation. 

In this paper we have only considered the production of the SM Higgs boson. 
The Minimal Supersymmetric extension of the Standard Model (MSSM) leads to 
two CP-even neutral Higgs bosons [1]. They are produced by gg fusion through 
loops of heavy quarks (top, bottom) and squarks. For small values (tan /3&5) of the 
MSSM parameter tan /3, the NLO QCD corrections to this production mechanism are 
comparable (to better than 10%) [12, 40, 9] to those for SM Higgs boson production. 
Therefore, the NNLO i^-factors computed in this paper could also be applicable to 
MSSM Higgs boson production. 

Note added: The calculation of the soft and virtual NNLO corrections to Higgs 
boson production has independently been performed in Ref.[41]. The method used 
in Ref.[41] is different from ours. The analytical results fully agree. 
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